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Abstract 

It is well-known that Birkhoff's theorem is no longer valid in theories with more than four dimensions. Thus, in 
these theories the effective four-dimensional picture allows the existence of different possible, non-Schwarzschild, 
scenarios for the description of the spacetime outside of a spherical star, contrary to four- dimensional general 
relativity. We investigate the exterior spacetime of a spherically symmetric star in the context of Kaluza-Klein 
gravity. We take a well-known family of static spherically symmetric solutions of the Einstein equations in an 
empty five-dimensional universe, and analyze possible stellar exteriors that are conformal to the metric induced 
on four- dimensional hypersurfaces orthogonal to the extra dimension. We show that all these non-Schwarzschild 
exteriors can continuously be matched with the interior of the star, indicating that the matching conditions at 
the boundary of a star do not require an unique exterior. Then, without making any assumptions about the 
interior solution, we prove the following statement: the condition that in the weak-field limit we recover the usual 
Newtonian physics singles out an unique exterior. This exterior is "similar" to Scharzschild vacuum in the sense 
that it has no effect on gravitational interactions. However, it is more realistic because instead of being absolutely 
empty, it is consistent with the existence of quantum zero-point fields (P.S. Wesson, Phys. Essays, Orion 5, 
591(1992)). We also examine the question of how would the deviation from the Schwarzschild vacuum exterior 
affect the parameters of a neutron star. In the context of a model star of uniform density, we show that the 
general relativity upper limit M/R < 4/9 is significantly increased as we go away from the Schwarzschild vacuum 
exterior. We find that, in principle, the compactness limit of a star can be larger than 1/2, without being a black 
hole. The generality of our approach is also discussed. 
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1 Introduction 



In Newtonian and relativistic theories of gravity the simplest models of isolated stars are provided by static spherically 
symmetric distributions of matter surrounded by empty space. In general relativity, the interior of a star is modeled 
by a solution of the Einstein field equations with some energy-momentum tensor satisfying physical conditions. The 
exterior of an isolated star is described by the Schwarzschild metric, which according to Birkhoff 's theorem is the 
unique asymptotically flat vacuum solution with spatial spherical symmetry. At the boundary of the star, which is 
a spherical surface of zero-pressure, both solutions have to be matched following a standard procedure in general 
relativity. These models are used to study the physical features of compact objects such as white dwarfs and neutron 
stars. 

Today, there are a number of theories that envision our four-dimensional world as embedded in a larger universe 
with more than four dimensions. In these theories it is natural to look for stellar models, similar to those of 
general relativity, and study the consequences of the existence of extra dimensions on the structure of neutron stars. 
However, there is an important aspect that complicates the construction of such models in these theories. Namely 
that Birkhoff 's theorem is no longer valid in more than four dimensions, i.e., there is no an unique asymptotically flat 
vacuum solution with spatial spherical symmetry. As a consequence, the effective theory in AD allows the existence 
of different possible, non- Schwarzschild, scenarios for the description of the spacetime outside of a spherical star, 
contrary to four-dimensional general relativity. 

In braneworld models, Germani and Maartens pQ have found two exact vacuum solutions, both asymptotically 
Schwarzschild, which can be used to represent the exterior of an uniform-density star on the brane. The analysis 
of the gravitational collapse and black hole formation on the brane is even more complicated; Bruni, Germani and 
Maartens [2] have shown that the vacuum exterior of a collapsing dust cloud cannot be static, although, on general 
physical grounds it is expected that the non-static behavior will be transient, so that the exterior tends to a static 
form. Similar results, for the collapse of a dust cloud, have been derived by Kofinas and Papantonopoulos [3] in the 
context of various braneworld models with curvature corrections. 

However, in Kaluza-Klein theories a milder version of Birkhoff's theorem is true. Namely, there is only one 
family of spherically symmetric exact solutions of the field equations Rab — which are asymptotically flat, static 
and independent of the "extra" coordinates. In five-dimensions, in the form given by Davidson and Owen |4J, it is 
described by the line element 

\ar + lj a 4 r 4 (ar - l) 2 k( fc - 1 )- 1 ] L 1 \ar - 1 J w 

where a is a constant with dimensions of L , and a as well as k are parameters that obey the constraint 

CT 2 (fc 2 -fc + l) = 1. (2) 

This family of solutions has been rediscovered in different forms by Kramer [S] and, although in a different context, 
by Gross and Perry [B] (for a recent discussion see [7] and references therein) . Its importance resides in the fact that 
in the limit k — > oo (ak — > 1), it yieldfQ 

?2_A^ 1 / flr \ j,2 ( i , 1 \ U „2 , 2 > n 2l i J„.2 



dS = h T~i — dt ~ U + — [dr 2 +r 2 dn 2 ]±dy 2 , (3) 

\1 + 1/ar J \ ar J 

which on every subspace y — constant reduces to the Schwarzschild metric, in isotropic coordinates, for a central 
mass M — 2/ a. Namely, 

In five-dimensional Kaluza-Klein theory, these solutions play a central role in the discussion of many important 
observational problems, which include the classical tests of relativity, as well as the geodesic precession of a gyroscope 



1 This metric is known as the Schwarzschild black string. 
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and possible departures from the equivalence principle [H] . In the context of the induced-matter approach, the metric 
is interpreted as describing extended spherical objects called solitons [SJ. In this interpretation the matter 
distribution contains a lightlike singularity at r = 1/a [7] [10j . which, in principle, can be visible to an external 
observer. Only in the black string limit k — > oo (ak — > 1) the hD metric ([1} possesses an event horizon; for any finite 
value of k the horizon is reduced to a singular point. 

However, the presence of naked singularities makes everybody uncomfortable. In order to avoid such singularities, 
we have to exclude the central region and require r > 1/a. What this suggests is that the asymptotically flat metric 
(JT]) should be interpreted as an "exterior" solution, describing the gravitational field outside of the core of a spherical 
matter distribution. In this interpretation, the effective exterior is not empty because there are nonlocal stresses 
induced from the Weyl curvature in 5D, which in AD behave like radiationo The "interior" region r < 1/a has to 
be described by some another solution of the field equations, which must be regular at the origin and not necessarily 
asymptotically flat. The aim of this work is to investigate this interpretation and its consequences on stellar models. 

In this context, the simplest candidate to describe the exterior of a spherical star is given by the spacetime part of 
the 5D vacuum solution (TTJ), which is the metric induced on every hypersurface y — constant. Nevertheless, we show 
here that this simple exterior leads to stellar models that have no Newtonian limit. On the other hand, according to 
the correspondence principle, we should require that our stellar models be consistent not only with general relativity 
but also with Newtonian models. Therefore, in this work we consider a more general scenario. The paper is organized 
as follows. 

In section 2 we use the five-dimensional line element ([T]) to generate a family of asymptotically flat metrics in 
AD which are conformal to the metric induced on y = constant hypersurfaccs. We examine this family as a possible 
choice for describing the four-dimensional spacetime outside of a spherical star of radius ri, > 1/a. 

In section 3, from the boundary conditions, we find that there is only one member of this family compatible with 
the weak-field limit. In other words, the requirement that in the weak-field limit we recover the usual Newtonian 
physics singles out an unique exterior. This exterior depends on one parameter, which we call e, and includes 
the vacuum Schwarzschild exterior for e — 1. The intriguing feature of this metric is that it is "similar" to the 
Scharzschild vacuum in the sense that it has no effect on gravitational interactions. However, it is more realistic 
because instead of being absolutely empty, it is consistent with the existence of quantum zero-point fields [llj . 

In section 4, in order to illustrate our results we examine in detail the particular model where the central core is 
represented by a perfect fluid sphere of homogeneous incompressible fluid. We calculate the physical parameters of 
a neutron star for distinct exteriors and different degrees of deviation from the Schwarzschild metric. We show that 
the general relativity compactness limit M/R < 4/9 is significantly increased as we go away from the Schwarzschild 
vacuum exterior. In fact, we find that in principle M/R can be larger than 1/2, without being a black hole. Finally, 
in section 5 we give a summary of our work and discuss the generality of our results. 

2 The model 

In this section we present the main steps of our work. First, we address the question of how to identify the effective 
four-dimensional spacetime. The correct answer to this question is crucial in order to be able to predict effects 
from an extra dimension. Then, from ([1]), we construct the "appropriate" set of possible exterior metrics. Next, we 
describe the model for the central core and examine the boundary conditions, which ensure that the interior and 
exterior metric join continuously. 

2.1 The effective spacetime 

In five-dimensional models, the effective four-dimensional pictur^l is usually recovered, or constructed, from gj^f 
the four-dimensional metric induced on AD hypersurfaces that are orthogonal to what is considered the "extra" 

2 It is sometimes called "black" or "Weyl" radiation, because in this case the induced energy-momentum tensor T^ v can be expressed 
as Tfiv = —tEfii/, where E^v represents the spacetime projection of the five-dimensional Weyl tensor, which is tracclcss. 
3 This is the picture registered by an observer who is not aware of the existence of an extra dimension. 
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dimension |12| . For a given line element in 5D 



dS 2 = g^{x",y)dx' l dx v + e$ 2 (x p , y)dy 2 , 



±1, 



(5) 



the induced metric, on any hypersurface y = constant, is just g^, i.e., g™ d = g^ v . The question is how to obtain 
our physical spacetime, with metric g e JJ , from the induced metric <7™ d . 

Unfortunately, there is no a consensus answer to this question. In fact, there are different competing approaches 
which are not equivalent. Today, in braneworld models, as well as in space-time-matter (STM) theory, the standard 
technique is to identify the metric of our spacetime with the induced one, i.e., g e JJ = g l ^f ■ However, there are a 
number of theories and approaches where g e JJ is constructed in a different ways. 

The simplest example is given by a metric which has been termed canonical [13j , where the g e JJ is conformal to 
9uv i w hh a conformal factor depending only on the extra coordinate. Also, a number of authors [8] [14] [15] [16] 
[T7] [T5] have considered models where g e JJ is conformal to g l ^f, with a conformal factor that is a "function" of the 
scalar field «3>. 

In this work we examine the consequences of AD models obtained from a dimensional reduction of the metric in 
5D, on a hypersurface y — constant, with a conformal rescaling 



[At/ 



(6) 



where N is an arbitrary real constant. This ansatz consolidates various approaches in the literature. In the context of 
Kaluza-Klein gravity and STM, it has been considered by Wesson [5j, Kokarev [H]-[TH], Sajko [IB] , and the present 
author [17] . For N = 0, it reproduces the usual interpretation of braneworld and STM theories. For N = —2 the 
interpretation is similar to the one provided by the canonical metric. For N = 1 it yields the classical Davidson-Owen 
[4] and Dolan-Duff [19] interpretation. 

2.2 Possible non-Schwarzschild stellar exteriors 

The effective four-dimensional spacetime (0, corresponding to the empty five-dimensional universe described by 
metric ([1]), can be written as 



ds 1 



where e = a(k - N/2) ancfl 



dt A 



(ar + l) 2 (ar-l) 2 



ar - 



2[e+<r(JV-l)] 



[dr 2 +r 2 dn 2 }, 



(7) 



e(l — N) + y/{N 2 - 2AT + 4) - 3e 2 



TV 2 — 2N + 4 



e(l — N) — y/(N 2 -2N + 4)- 3e 2 



N 2 - 2N + 4 



for N > 1: 



for N < 1. 



(8) 



This parameterization is useful because the "deviation" from the exterior Schwarzschild metric can be measured 
by |e — 1|. Since the exterior spacetime is not in general empty, from a physical point of view taking e 1 is 
equivalent to giving the star an "atmosphere" extending all the way to infinity with a matter density decreasing as 
p ~ l/a 2 r 4 . The special case where the star has a sharp surface corresponds to e = 1, where we recover the exterior 
Schwarzschild metric, for any value of N. In Appendix A we present the equation of state and gravitational mass, 
of such an atmosphere. Thus, the questions of interest here are the following: 

1. How does the existence of stars restrict the possible values of TV? 



4 We note that for N = 1 the choice of a is immaterial, because the metric depends only on e. However, we have selected a positive a 
in l(8} for consistency with our previous work [7]. 
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2. How does a possible deviation from the Schwarzschild vacuum exterior, i.e., going away from e = 1, can affect 
the star parameters? 

The total gravitational mass, measured at infinity, is given by 

M = -. (9) 

Thus, in the Schwarzschild limit M = 2 /a. Therefore, to ensure the positivity of M, we will choose a > and e > 
everywhere. 

We note that a(N — 1) is invariant under the transformation N — > (2 — N). This means that the effective 
spacetime for N = is identical to the one for N — 2. The same for N = 3 and iV = — 1; N = 4 and N = —2; and so 
In other words: the effective spacetimes with N > I duplicate the ones with N < 1 . Only the one with N = 1 
remains the same. 



2.3 The interior solution and matching conditions 

In this work we restrict our discussion to the case where the interior solution is time independent. The matching 
conditions are easier to derive if we take the line element inside the fluid sphere with the same symmetry properties 
as the exterior metric ([7]). Namely, the interior static will be taken as 

ds 2 = e v{r) dt 2 ~ e x{r) [dr 2 + r 2 {d9 2 + sin 2 6d(j) 2 )]. (10) 

The three-dimensional surface S, defined by the equation 

r-r b = 0, (II) 

where is a constant, separates the spacetime into two regions: r > and r < r^, described by the metrics ([7]) and 
(fTU)l . respectively. 

The problem of joining two distinct metrics across surfaces of discontinuity has been discussed by a number of 
authors (see [5D] and references therein). Robson [5D] showed the complete equivalence of the junction conditions 
proposed by Lichnerowicz and by O'Brien and Synge at non-null hypersurfaccs. The latter are analogous to the 
requirement that the first and second fundamental forms be continuous at S. 

For the line elements under consideration, the first fundamental form at E, becomes 

ds h--r=r b ) = e u{rb) dt 2 - e x(r ^r 2 (d0 2 + sin 2 9 d(f> 2 )}. (12) 
Therefore, its continuity across E requires 

e u(r b ) = f ar b -l \ 2£ 
\ar b + l) ' 

eA(? . b) = (ar b + l) 2 (ar b ~l) 2 Ur b + 1 \ 
a 4 r^ \ar b -l) 

Besides, the "physical" radius R b of the sphere with coordinate radius r b is given by 



(13) 



A =,, e ^=(^(^y• + ' < ~- l,, . (14) 



a 2 r b \ar b — 1 
The second fundamental form, say dK 2 , at E is given by 

dK 2 , = n^dx^dx", (15) 



5 These are just examples, we are not assuming that N is an integer number 
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where n M is the unit vector orthogonal to S. In the present case, if we take 



then we find 



n„ = c^ 2 , (16) 



_ _ -A(n)/2 



2 V / Vr,, 



(17) 



Consequently, in isotropic coordinates, the requirement of continuity of the second fundamental form across £ is 
equivalent to demand that the first derivatives of the metric be continuous across^ S. 
For the exterior metric ([7]), the continuity of v' and A' yields 

v (r b ) = 



(a^rf-l)' 

4{1- [g + tr(JV-l) ] ar b } 
r b (a 2 rl - 1) 



A (n) = Z-7Z2T2 7^ • ( 18 ) 



From a physical point of view, the continuity of v' is equivalent to the continuity of the gravitational mass (|A-11I) , 
in such a way that the force on test particles located at the surface will be the same whether we calculate it from 
inside or outside. While, the continuity of A' means that the radial pressure (i.e. — 2f) of the fluid at the surface 
joins continuously with the pressure of space outside the sphere. 

Instead of (|10|) . we may consider the most general static interior metric, which in curvature coordinates can be 
written as 

ds 2 = e^dt 2 - e^dR 2 - R 2 (d8 2 + sin 2 Odtf ). (19) 

Now the continuity of the second fundamental form does not require the continuity of dfi/dR. However, by virtue of 
the equivalence between Lichnerowicz's and O'Brien-Synge's junction conditions, this and the exterior metric match 
if the mass and radial pressure are continuous at E. 



Embedding the interior solution in 5D: The interior solution (fTQ| . as well as (|T9|) . can be embedded in the 
five-dimensional world by foliating the 5-dimensional manifold as a set of 4-dimensional hypersurfaces (or slices), in 
such a way that the effective metric in AD coincides with the induced one on every y — constant hypersurface. The 
embedding can be taken as 

dS 2 = e v{r) dt 2 - e A(r) [dr 2 + r 2 (d6 2 + sin 2 0d(/) 2 )} ± dy 2 . (20) 

which physically corresponds to a star carrying no-scalar charge [3T] , [32] . The corresponding field equations in 5D 
are 

G AB = k 2 i5) ^T AB , (21) 

where k 2 ^ is a constant introduced for dimensional considerations, ab is the energy-momentum tensor in 51? 
and the non- vanishing components of the Einstein tensor Gab are 



Gn = e 



G\ = -e 



A 



A" + 


A' 2 
T + 


4). 


A' 2 


AV 










A" + 


v" 

h 


v' 2 A' + v' 






4 + 2r 



Gi = -e~* l -^ + ^r + ^— , (22) 



3 Thus, in isotropic coordinates the "weaker" junction conditions of O'Brien and Synge are equivalent to those of Lichnerowicz. 
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G 3 3 =Gl and 

G\ = -e- x A" + — + 2— + — + — + — + — . 23 

4 V 4 r 4 r 4 2 / v/ 

We note the relationship Gg + G\ + 2G\ — 2G|. Thus, taking ^T B — diag(p, — p r , —p±, —p_i_, —Pb) (see equation 
(|A-5p ) we find that the matter in 5D satisfies the equation of stat^ll 

p = p r + 2rl - 2p 5 , (24) 

which is the counterpart of the equation of state (|A-8j) for the exterior of a spherical star. A four-dimensional 
observer living in a subspace y = constant is not directly aware of (|23p. the "pressure" along the extra dimension, 
and therefore will not be able to find such a simple equation of state between the density and pressures. However, 
this does not mean that she/he loses all information about the 5th dimension; it is consolidated in the nonlocal 
stresses induced in 4D from the Weyl curvature in 51?. As a consequence, even in the absence of matter, the exterior 
of a spherical star is not in general an empty Schwarzschild spacetime. We note that in the above discussion the 
signature of the extra dimension is irrelevant. 

An observer in 4D interprets the spacetime part of ^Tab as the energy- momentum in 4D. More precisely, as 

8irT a/ 3 — G a p. (25) 

It is clear that the five-dimensional conservation equation ^T^-b = 0> m the present case yields the conserva- 
tion equation = in 4D, which is equivalent to the "generalized" Tolman-Oppenheimer-Volkov equation of 
hydrostatic equilibrium 

(P + Pr)^e-^V> = ~f+(l + X') (P±-Pr). (26) 



r 



Finally, we note that for metrics (flT))) and (|2T)|) a motion along a geodesic in AD is also geodesic in 5D. Therefore, a 
particle located on a hypersurface y — constant, with zero momentum along the extra dimension, will stay on that 
hypersurface. 



3 General interior and the weak-field limit 



The main question in this paper is whether the existence of stars leads to some constraints on N. In this section we 
give a positive answer to this question. We use the boundary conditions, as well as the weak-field limit, to show that 
the exterior of a spherical star is given by the metric jJJ with N = 1; for any other N ^ 1 there are problems in the 
Newtonian limit. 

First of all, let us introduce the notation 

x = ar^, 



(27) 



in terms of which the continuity equations (|18p become 

v'{r b ) _ e 



4a 
A'(r fe ) 
4a 



x > 1, 



x 2 - 1 

{1- [e + <7(JV-l)] x} 
x (x 1 — 1) 



These can be used to get an equation for x. Namely, 



A'(r b ) 



V{r h ) 



-ct(I-JV) 



(28) 



(29) 



7 Here p,p r and pj_ represent the energy density, radial and tangential "pressure", respectively; ps is the pressure along the extra 
dimension. 
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In the limit s — > 1, the radial pressure vanishes at the boundary surface, i.e., we recover the usual matching conditions 
for the Schwarzschild exterior. However, we note that in the weak-field limit (|2"5)) violates the positivity of x. Indeed, 
in the Newtonian limit (see equation ([55)1 bellow) 

£—1. (30) 

On the other hand <r(l — N) > for any N ^ 1 and £ ^ 1. Therefore, X'/v' cannot approximate —1 otherwise x < 0. 
This shows that, given the standard matching conditions at the boundary, the effective exterior of a spherical star 
must be described either by the vacuum Schwarzschild spacetime or by the line element ([7]) with N = 1, otherwise 
the model star is inconsistent with a Newtonian limit. 

Newtonian limit: In the weak field- limit we set e v = 1 + £ f(r) and e A = 1 + £ h(r), where £ is a "small" 
dimensionless parameter, i.e., |£| <§; 1. In this limit p^$> p. Therefore, from (|A-10[) it follows that 

M(r) = 4tt / f 2 p(f)df + 0{£ 2 ). (31) 
Jo 

In this limit the energy density is given by 

8irp = -£ (h" + ^pj . (32) 
From which we get M — — £r 2 /i'/2. On the other hand, in this limit, (|A-11|) yields M = £r 2 /'/2. Consequently, 

f'(r) = -h'(r), - = -! + £(/!-/). (33) 
Then, the pressures are of second order in £. Namely, 

8Trp r = (8h'h + rti 2 ) + ? (f + h), 

f 2 f 2 (rh" + h') 

8np ± = -1- (3rti 2 + 4rfrfr" + Ahh') + ^ 1 " ' (/ + fe). (34) 

We note that 



87T 



dp r 2 

-r- + -(Pj. "Pr) 
ar r 



^ + (35) 



The r.h.s. can be written as — (£h' /2)(8irp) — 8irpM/r 2 . Thus, we recover the Newtonian equation of hydrostatic 
equilibrium, viz., 

Mp _ dp r + 2 ( 
r 2 c?r r 

which is nothing but the non-relativistic limit of the TOV equation (f26|) 



■^2-=-^T+-(P-J--Pr)' ( 36 ) 



The non-Schwarzschild stellar exterior: The conclusion from this section is that, in Kaluza-Klein gravity the 
only exterior spacetime that is consistent with the existence of a Newtonian limit is the one given by jJJ with N = 1. 
In this case, the exterior of a spherical star takes a particular simple form. Namely, 

In addition to the equation of state p = —{p r + 2pj_) mentioned in (|A-9[) . we find that the external field obeys a 
second equation of state, which is 

P± = ~Pr- (38) 
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In other words; p r = p and p± = — p, where the energy density is given by 



4«V(l-6 2 ) (ar-lV* 
8lTP= {ar + mar-iy {^TT ) ' (39) 

and the gravitational mass (|A-12|1 is constant outside the fluid sphere. Namely, 

M =—. (40) 
a 

This means that the (total) gravitational mass resides entirely inside the fluid sphere; the external cloud has no 
influence on the gravitational mass. This is due to the equation of state p — — (p r + 2pj_), which describes matter 
with no effect on gravitational interactional Finally, we note that the conditions p > and M > demand < e < 1. 
For e = 1 we recover the Schwarzschild exterior. 



4 Model for the stellar interior 

In order to illustrate the above results we need to assume some interior model. Here we consider a model star in 
the shape of a homogeneous sphere of incompressible fluid with radius r&. Thus, inside the sphere the metric will be 
the interior Schwarzschild metric, and outside ri, it will join continuously to the exterior metric ((7|). Despite their 
simplicity, uniform-density models give insight about much wider classes of stellar models |26] . 

The line element describing the interior of a fluid sphere with constant energy density po and isotropic pressures 
can be written as 

The matter quantities are 

12C 

8ir Po = —, <7>0, (42) 

8»p« = — mi+m — • (43) 

In order to simplify these expressions we evaluate the pressure and the equation of state at the center. Namely, 

4(5 -2(7) 

8nPc= D2 -, Pc = lPo, (44) 

where p c = p(0) and 7 is a constant. In terms of these quantities we find: B = C(2 + 37); (B — 2C) = 3Cj and 
(C — 2B) = — 3C(1 + 27). Consequently, the appropriate line element and pressure are 

dg2 = A2 / l + (2 + 3 7 )(7r 2 \ 2 a _ 3C 2 ^ a 

V 1 + Cr 2 / 27rpo(l + Cr 2 ) 2V ' v ' 

[ 7 -(l + 2 7 )Cr 2 ] 

P = / > + (2 + 37)Cr 2 ] - (46) 
Since the density and pressure are positive, it follows that 

7 > 0. (47) 

We note that, in principle, 7 can take arbitrarily small positive values; there is no other restriction. 



8 This equation of state is the generalization of p = — 3p, which has been considered in different contexts by several authors. Notably, 
in discussions of cosmic strings by Gott and Rees |23| . by Kolb 1241 . as well as in certain sources (called "limiting configurations") for 
the Reissner-Norstrom field by the present author I2r>| . Besides, it is the only equation of state consistent with the existence of quantum 
zero-point fields by Wesson 
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For completeness, we note that although not directly observed in AD, the pressure along the extra dimension p§, 
can be obtained from the equation of state (|2"4")l as 

p 5 = \(3p-p). (48) 

Then from (|46|) we find 

[(3 7 -l)-(5 + 9 7 )Cr 2 ] 

P5=P0 2[l + (2 + 3 7 )CV 2 ] ■ (49) 

4.1 Boundary conditions for the Schwarzschild interior 

In terms of the dimensionless quantities £, x and Y defined as 

C = C/a 2 , x = ar b , Y — aRb, (50) 

we have Cr\ = (x 2 and 

With this notation the continuity of the first fundamental form gives 



[1 + (2 + 3 7 )Ca; 2 
and 



x z - 1 



xY 



(52) 



From the continuity of the second fundamental form we get 

(3 7 + 1)C* 



Y= /n , m ^ (53) 



(l + Ca; 2 )[l + (2 + 3 7 )C2: 2 ] (x 2 - 1) ' 

and 

(x 2 _[e + a(N - l)]x - 1 



1 + Qx 2 {x 2 - 1) 

These equations require x > 1 and [e + cr(./V — l)]a; ^ 1, otherwise C = 0, for all values of e and N. 

4.2 Solution of the boundary conditions 

Dividing (f54"| by (f5"5"l) we obtain 

3 7 + 1 ex 



(54) 
(55) 



[l + (2 + 3 7 )(Cx 2 )] [e + cr(iV- l)]x- 1' 

from which we get 



(56) 



, [(iV — 1)(3 7 + 1)ct + 3e 7 ] a? - (1 + 3 7 ) 



On the other hand, from ([5 



^ 2 )=^ — v £ ; (2+ 3 7 i J v — - w 

2 [ g + f 7(jV-l)]x-l 
(C * } = x[x-e-a(N-l)} - (58) 



We note that x ^ e + cr(A r — 1), otherwise for e = 1 one would have a; = 1. 
Equating (|57|) and ((58)) we get the desired equation for x. Namely, 

[(l + 3 7 )(7V-l)cr + 3e 7 ] x 2 -(l + 3 7 )[(iV-l) 2 ( T 2 + 3(iV-l)ecr + 2e 2 + l] x+ [(1 + 3 7 )(7V- l)cr + 3e(l + 2 7 )] = 0. (59) 
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Physical conditions: The solutions of the above equations have to satisfy physical requirements. Positivity of 
mass requires e > 0; positivity of po requires C > 0, which translates into £ > 0. Consequently, from (|57p we get 

Po > => CI := [(iV - 1)(3 7 + 1)<7 + 3e 7 ] > 0, (60) 

which also implies that the positivity of the coefficient in front of x 2 . On the other hand, from (|A-1[) it follows that 

p > =>■ C2 := [37V(2 - N)a 2 + 4e(l - AT) a] > 0. (61) 

4.2.1 Solution for e = 1 

In order to facilitate the discussion, we first discuss the solution of the boundary conditions for the Schwarzschild 
exterior metric. This will allow us to compare and contrast the parameters of our model star obtained under different 
conditions. 

In this case a = 0, and the above equation for x reduces to 

jx 2 - (1 + 3 7 )x + (1 + 2 7 ) = 0, (62) 

whose solution is 



l + 2 7 
x = . 

7 

The other "solution" x = 1 is excluded. Substituting this expression into (|58[) we find 



(63) 



and from (|5ip we obtain 



(C 2 ) = ^ (64) 



r = fet£. (65) 



7(1 + 27) 

Now, the original constants can be obtained from ([50]) and ([52")) . Namely, 

1 3C /l + 3 7 \ 6 ^ 2 7 (r^ 1 



In addition, we find 



A = = , Cr 2 = '— - . (66) 

1 + 37' 2^0 \i + 2 7 y i + 2 7 v^y v ; 

A/(l + 2 7 ) /o«rn~M 2V3 7 3 /2(i + 2 7 )^ 

where the total mass has been obtained from (|53[) - Thus, the boundary conditions allow us to obtain the value of all 
the constants in the solution in terms of the equation of the state at the center. 



Buchdahl limit: In the present case, from ([9]) it follows that a = 2/M. Consequently, using (|65|) we obtain 

M _ 27(1 + 2 7 ) 
W b ~ (3 7 +l) 2 ' 

In the limit 7 — > 00 we get the famous absolute upper limit 



(68) 



^f) = I « 0.444, (69) 

K bj max y 

discovered by Buchdahl [27], for all static fluid spheres whose (i) energy density does not increase outward and (ii) 
the material of the sphere is locally isotropic. 
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Newtonian limit: In the limit 7^0, from (|66p we see that A — > 1 and (3C/2irpo) — » 1. Also, since < (r/rb) < 1, 
it follows that Cr — > 0. Consequently, e v — > 1 and e A — > 1. Therefore, in this limit (|2l)|) reduces to 



which is the usual balance equation in Newtonian physics. 

4.2.2 General solution: e ^ 1 

The general solution of (|5"9")) can be written as 

[(N - l) 2 fi 2 + 3ea(N - 1) + (1 + 2e 2 )] + y/F 



(70) 



2(iV- l)cr + 67e/(l + 37) 



(71) 



where 

F = (N — 1)V + 6(A^ - l) 3 ea 3 + (TV - l) 2 (13e 2 - 2)<t 2 + 6{N - l)(2e 2 - l)ea + (2e 2 - l) 2 + 4 ± . (72) 

The sign of the root in (|71[) has been chosen positive to make sure that in the limit e — > 1 we recover the solution 
discussed in the preceding subsection. 

Thus, for any given (iV, e, 7) we get cc from (f7T|) - Substituting the result into ([55f we find (C^ 2 )- Using them in 
(1ST))) we obtain the constants C, r;, and A. Equation (f6"6"|) epitomizes the fact that the first two constants depend on 
the total mass. In the present case the total mass of the configuration, i.e., the mass of the fluid sphere added to the 
mass of the external cloud, is obtained from © and (|53")l . The result is 

2 V /[3(C^ 2 )] e x(x - l)U e+,T ( N -V-y 



y^p° M = ( i+c^)( a: +i)[(^(iv-i) + i] ■ ( 73 ) 

Warning M 7^ eMschw- At this point, in order to avoid a misunderstanding, the following comment is in order. 
From ([9]), one could conclude that, since for e = 1 (which corresponds to the Schwarzschild exterior metric) the total 
mass is Ms c hw = 2/a, then the total mass for any e ^ 1 is just M = eM$ c hw However, this is not true; it is easy 
to see that (|73M ex ([57]) . This is because a is not a fixed constant but it is given by the boundary condition (|53[) . 
What is true is that for e = 1, and the corresponding values of x and (C,x 2 ) given by ([63]) and ([64]) . we recover the 
mass in the Schwarzschild case. 

Coming back to our discussion we calculate the mass of the star, say Mb. By virtue of the continuity of the 
second fundamental form, it can be easily obtained from (|A-12I) . Namely, 

lX o-(l-iV) 



Mi = M I — ^— j- J . (74) 



We note that Mb = M only in two cases: for the Schwarzschild exterior metric (a = 0), and for N = 1. 
Now, the gravitational potential at the surface of the star is given by 

M b _ 2ex(x - l)^- 1 ) 

£^ ~ (lE + l)^ 1 ) 



(75) 



For practical reasons, it is useful to express the mass of the star in terms of the solar mass M@, and its radius in 
kilometers (km). They are given by 



1.089 x 10 9 JlCx 1 ) e x(x - l)^" 1 ) 
M b = == — — - i - x M , 

^/R(i + i)^)[i + (Cx 2 )] 

0.803 x 10 9 7(C^) 
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where [po] denotes the numerical value of the energy density in (gr/cm 3 ). 

For e = 1 the star is surrounded by empty space with p — p r = p±_ = for r > r\>. However, for any s ^ 1 it 
does have an atmosphere, whose matter density decreases as p ~ l/a 2 r 4 from its maximum value, p\s+, which is just 
outside S : r = rfc, to zero at infinity. For the case under consideration, from (|53p . (|A-1[) and using (158p we find 

= xa 2 [3N(2~N) + 4e(l~N)/a] 

po 12 [x - e - a(N - l)][x{e + a(N - 1)) - 1]' 1 ' 

In view of the complexity of the general solution, we will analyze it for different values of e and ./V separately. 

4.3 Factorization with N = 1 

In this case the above expressions are specially simple. As in the Schwarzschild case, we find that for any fixed £ / 1 
the radius and the mass of the fluid sphere are very sensitive to the equation of state 7. This is because x decreases 
very fast, as x ~ 1/7, for small 7, to a constant value for relatively large values of 7. 

4.3.1 Limiting configurations 

Let us first examine the solution of (fTTj) in the limit 7 — > 00. In this case from (j72|) we get F = (2e 2 — l) 2 . Since 
we are looking for solutions that are continuously connected to the Schwarzschild one, we take e 2 > 1/2. With this 
choice, from (I7T|) we obtain 

x = 2s. (78) 

Then from we find 

2e 2 - 1 

(C^H^-A (79) 

For e = 1 we recover x = 2 and ((s 2 ) = 1/2, which are the corresponding limiting values in the Schwarzschild case. 
Using ([78]) in ([S0 ]) -([S3 f we find the constants 

r _2s 2 -l A . (2e-l)(^) 

2e 2 M 21 3 7 (2e 2 -l)(2£ + l)(^ 1 )' 1 j 

as well as the mass M, as a function of e. Namely, 



M 4 £ y[3(2e 2 -l)A] (2 £ -l)(^ 2 ) 

M = v» + D (e+2) ' (81) 



From ([5T|l and ([40]) we get the upper limit of the surface gravitational potential 



Buchdahl limit for different e: For e = (1.00; 0.80; 0.75) we find 

VpoAf = (0.145; 0.135, 0.112), 

~p^R b = (0.326, 0.265; 0.207), 
M 



- 



Rb 



4/9; 0.508; 0.538), (83) 



respectively. Thus, with the decrease of e both the mass and the radius of the star decrease, in such a way that its 
surface gravitational potential increases. We note that M/Rb > 1/2, for certain values of e. 
The red-shift Z(r) of the light emitted from a point r in the sphere to infinity is given by 



Z(r) = l/ v / toW-l. (84) 



13 



In the present case, the red-shift of the light emitted from the boundary surface is given by 



Since, in the present lining], e 2 > 1/2 it follows that the maximum surface red-shift is Zb = 2.478. This indicates 
that, although the surface gravitational potential is M/Rb > 1/2 (instead of the usual M/Rb < 1/2), the star is not 
a black hole. 



Newtonian limit for different e: For 7^0, from (|7Tj) with N — 1, we get 

2e 2 + l 

Substituting this into (|5"T)) . or we find 



(86) 



For £ = 1 these quantities become x — * 1/7, (C^ 2 ) * 7> which are identical to their Newtonian values in the 
Schwarzschild case. In this limit, (fSTj) yields Y ~ x Substituting these expressions into (f52"| and (j53"|) we get 

A-l, — -1, Cr 2 ->4^_f^V (88) 
2^ 2e 2 + 1 V^/ 

Therefore, for ./V = 1 in the limit 7 — > we recover flat spacetime as well as the Newtonian balance equation (|7D|) . 



4.3.2 Effects of e 7^ 1 on neutron star parameters 

The second important question in this work is how going away from e = 1 may affect the star parameters. In order 
to study this question, we have calculated the surface gravitational potential M/Rb, the red-shifts at the center Z c 
and at the boundary Zb for different values of 7. The mass and radius of the fluid sphere are given by (|76]) . To 
obtain some numerical values for M and Rb we have taken p = 2 x 10 1A gr /cm? which has been used in different 
models of neutron stars [5S], [5^. Table 1 shows the results for the Schwarzschild solution (e = 1), which we use to 
compare and contrast with the results of Tables 2 and 3 for e = 0.9 and e = 0.8, respectively. 



Table 1: Star parameters for the Schwarzschild exterior 


7 =Pc/Pc 


M/Rb 


z c 


z b 


Rb(km) 


M g (M G ) 


1/10 


0.14 


0.30 


0.18 


15.13 


1.46 


1/3 


0.28 


1 


0.50 


21.16 


3.99 


1 


0.38 


3.00 


1.00 


24.59 


6.25 


3 


0.42 


9.00 


1.50 


26.02 


7.41 


10 


0.44 


30 


1.82 


26.54 


7.87 


00 


4/9 


00 


2.00 


26.77 


8.08 



Table 1 shows how the physical features of a neutron star, with a vacuum Schwarzschild exterior, depend on the 
equation of state at the center. 

9 In order to avoid misunderstanding, it should be noted that for finite values of 7, like 7 = 1/3 or 7 = 1, the parameter e is not 
restricted by this condition. 
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Table 2: Star parameters for e = 0.9 and N = 1 



7 = .Pc/Pc 


M/Rb 


Z c 


z b 


Rb{km) 


M b (M Q ) 


P|£+/P0 


1/10 


0.13 


0.28 


0.17 


14.56 


1.32 


0.007 


1/3 


0.28 


0.91 


0.47 


20.34 


3.71 


0.024 


1 


0.37 


2.68 


0.95 


23.57 


5.99 


0.061 


3 


0.43 


7.88 


1.48 


24.85 


7.27 


0.117 


10 


0.46 


25.88 


1.85 


25.26 


7.82 


0.167 


oo 


0.47 


oo 


2.09 


25.40 


8.07 


0.204 




Table 3: Star parameters for e = 0.8 and N = 1 


7 = Pc/Pc 


M/R b 


Z c 


z b 


Rb(km) 


M b (M ) 


P|E+/P0 


1/10 


0.13 


0.28 


0.17 


14.56 


1.32 


0.007 


1/3 


0.26 


0.80 


0.42 


19.25 


3.34 


0.054 


1 


0.37 


2.26 


0.87 


22.04 


5.54 


0.149 


3 


0.44 


6.22 


1.41 


22.71 


6.82 


0.333 


10 


0.48 


19.62 


1.86 


22.71 


7.34 


0.576 


oo 


0.51 


oo 


2.23 


21.79 


7.50 


0.857 



Tables 2 and 3 contain the discontinuity of the energy density at the boundary calculated from (|77|) . We see that as 
the surface gravitational potential increases, also increases, that is the atmosphere gets more concentrated near 
the surface, which is physically reasonable. 



4.4 Factorization with N ^ 1 and e ^ 1 

The energy condition (|60p requires the denominator in (|7ip to be positive. However, it vanishes at 7 = 7 given by 

(1 - N)a 



7 



3[e+(l-JV)<7] 



(89) 



We note that 7 = only for TV = 1 and the Schwarzschild exterior. Since a is positive for N < 1, and negative for 
N > 1, it follows that (1 — N)a > 0. Consequently, for any other N and e we have 7 > 0. 

On the other hand, since x is positive for 7 > 7 and negative for 7 < 7, it follows that 7 has a positive lower 
bound, i.e., 

7 > 7. (90) 

Thus, in accordance with our discussion in section 3, we conclude that the factorization with N 7^ 1, and e 7^ 1, 
is incompatible with the Newtonian limit, which requires 7 — > 0. Let us consider some particular cases with more 
detail. 



4.4.1 Factorization with N = 

This corresponds to the case where the metric of the physical spacetime is identified with the metric induced in 4D, 
which is the typical approach in induced-matter and braneworld theories. In this case the exterior metric ([7]) can be 
written as 

*. . r^yv - <"- +i >:'r" 2 {"-L±±t^ ] [*»+*«"], (9i) 

\ar + 1 / a 4 r 4 \ar — 1/ 

where we have used that, from ((SJ), in the present case 



e - V4 - 3e 2 

a = . 

2 

which requires e < 2/V3- 



(92) 
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Now, from ([50)1 and (fBTj) we get the constraints 



CI := -(37+ l)er + 3£7 > 0, 

C2 := 4ecr > 0. (93) 
From which we find that the positivity of the energy density inside of the fluid sphere requires 



e - V4 - 3e 2 

^= 3(« + v4-^) - ,94) 

We note that x is positive for 7 > 7 and negative for 7 < 7. In order to get a better understanding of these 
constraints, let us take e = (0.8, 0.9, 1, 1.1, 1.15, 2/^). We find 

7 > (-0.09, -0.05, 0, 0.09, 0.23, 1/3), 

a = (-0.32, -0.18, 0, 0.25, 0.48, 1/V3), (95) 

respectively. Thus, the positivity of the energy density outside of the fluid sphere requires e > 1. For e = 1 we 
recover the usual Schwarzschild picture, as expected. However, for e > 1, 7 has a lower positive bound, e.g., 7 = 0.23 
for e = 1.15, which increases with the increase of e. Certainly these models have no Newtonian limit. 

4.4.2 Factorization with N = 2 

In this case 



- + ^f!. (96) 



2 

Substituting this into ([7]) we find that the line element is identical to (|9ip . which is consistent with the symmetry 
under the change N — ► (2 — iV) noted in subsection 2.2. Also, both factorizations, A = and N = 2, yield the same 
equation of state (|A-6[) . and have the same lower limit on 7. 

The discussion in this section epitomizes a situation frequently found in physics. Namely, that the analysis of 
some particular "simple" example can be much more involved, in technical details, than the general case. However, 
the effort is worthwhile: the main conclusions from our uniform-density model are consistent with the general results 
presented in section 3. Specifically, (i) Only the factorization with N = 1 produces stars models that have a 
Newtonian limit; (ii) the factorizations with A and (2 — A) are physically indistinguishable. 



5 Summary and concluding remarks 

For theories in more than four-dimensions the crucial question is, how to recover the physics of our 4-dimensional 
world. The answer to this question is important in order to be able to calculate and predict specific effects from 
extra dimensions. For example, the possible existence of extra "non-gravitational" forces on test particles and their 
effect on the classical tests of relativity and the principle of equivalence [8] , [30] . 

In this paper, we have considered this question from the point of view of the stellar structure. Namely, since 
Birkhoff 's theorem is no longer valid in more than four-dimensions, there is no an unique asymptotically flat vacuum 
solution with spatial spherical symmetry. In other words, regardless of whether we choose to work in Kaluza-Klein 
gravity [9], STM [5], braneworld models pQ, [2] or others theories with curvature corrections [3], the fact is that 
the effective four-dimensional theory allows the existence of different possible, non-Schwarzschild, scenarios for the 
description of the spacetime outside of a spherical star, contrary to four-dimensional general relativity. 

In five-dimensional Kaluza-Klein gravity the well-known solutions due to Kramer-Davidson-Owen-Gross-Pcrry 
seem to represent the natural generalization of the Schwarzschild spacetime |9]-[12]. Even in this simple case of 
spherical symmetry, in ordinary three space, there are various possible options leading to different exteriors for a 
spherical star. A popular assumption is that we recover our AD spacetime by going onto a hypersurface y = constant. 
With this assumption, the information about the fifth dimension is consolidated in the nonlocal stresses induced in 
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AD from the Weyl curvature in 5D. As a consequence, even in the absence of matter, the exterior of a spherical star, 
like our sun, is not in general an empty Schwarzschild spacetime. 

However, early in our research we found that this simple approach leads to contradictions in the Newtonian limit. 
Therefore, in section 2 we generalized our discussion by considering, as a possible exterior for a static spherical star, 
the family of asymptotically flat metrics represented by ([7]). These are conformal to the induced metric ([6]) and 
constitute a generalization to a number of approaches in the literature (see for example 13J- [1 9 ) . 

The main question under scrutiny here was whether the existence of stars can be used to constraint the number 
of possible candidates for exteriors. Using the junctions conditions at the boundary of the star, and without making 
any assumptions about the interior metric, we have showed that (|37[) is the only exterior metric compatible with 
both; the weak-field Newtonian limit as well as the general relativistic Schwarzschild limit. All the other metrics 
with N yt= 1 do not have a Newtonian counterpart, although they are consistent with the Schwarzschild limit. 

It should be noted that the resulting effective matter (|3"9")) . outside of a spherical star, is "similar" to Schwarzschild 
vacuum in the sense that it has no effect on gravitational interactions. In fact, it satisfies the equation of state 
p = -~(p r + 2pj_), which as we can see from (|A-10[) corresponds to M = 0. This is the generalized version of the well- 
known equation of state p = —3p for isotropic fluid, which has been discussed in different contexts (see for example 
[23]-[ll]). One can speculate that this exterior (with N = 1) is more realistic than the Schwarzschild one, because 
instead of being absolutely empty (p = p r = p± = 0), it is consistent with the existence of quantum zero-point fields 

Em- 

The question may arise about the generality of our result. Namely, what would happen if we assume that the 
effective four-dimensional metric is like 

g$t = F(<%™ d , (97) 

where F is some "function" of $?. It is not difficult to show, using the boundary conditions, that the only function 
compatible with a Newtonian limit is F ~ $, i.e., the factorization with N = 1. 

Another important question here is what are the consequences of going away from the Schwarzschild exterior 
metric. Using a uniform-density model for the interior of the star, in Tables 1-3 we show how some physical quantities 
as the mass, radius and surface gravitational potential of a star are affected by an extra dimension. For the same 
equation of state at the center, we see that both the mass and the radius of a the neutron star are reduced as we 
separate from the Schwarzschild vacuum exterior. The most striking feature is that, in principle, as a consequence 
of the extra dimension, stars may exist with gravitational potentials larger than 1/2 without being black holes. 
Certainly, this can provide some clues for the observational test of the theory. 

Regarding the classical tests of general relativity and other experimental observations, within the context of 
Kaluza-Klein gravity [8], they have been calculated under the assumption that the physical spacetime outside an 
astrophysical object, like our sun, was described by the spacetime part of {2}. Our results here suggest to consider 
the metric (|3?]l 

ds * = 2e & . J_ (g + 1 ) 2[£+1] ^ + r wi 

As far as we know this line element has never been used in this kind of discussions. 



Appendix A: Possible non-vacuum exteriors of a spherical star. 

An observer in 4D, who is not directly aware of the existence of an extra dimension, will interpret the metric (O as 
if it were governed by an effective energy-momentum tensor. 

= aW[ZN{2-N)+te(l-N)/a] Ur-l \ 2 ^ N '^ ^ 



(ar + l) A (ar-l) 4 \ar + 1 

, Aa 5 r 3 a 2 Ul~N)[ar{Aa-2e-Na)/(2a) + {a 2 r 2 + l)/a] + {3/A)N 2 ar} far-l\ 2l£+a{N ' 1)] , , 
8 " Tl= (ar + i)4 (ar _ 1)4 (^T+iJ > ( A "2) 
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2 2a 5 r 3 a 2 {2ar[(2e + Na)(l - N)/(2a) - 1] - (1 - N)(a 2 r 2 + I) /a) + (1/2) A^V} / ar _ 2[ e + CT (Ar-i)] 
8TrT 9 = 1 



(ar + l) 4 (ar - l) 4 V fflr + 1 

(A-3) 

and Tg = Tf , as expected by virtue of the spherical symmetry. These equations show that T 2 2 7^ f° r an U value 
of AT and e ^ 1. Thus, the effective matter behaves as an anisotropic fluid, which can be described by an effective 
energy-momentum tensor of the form 

T^v = (p + P±)u^u v - P±9^ + (Pr - P±)x^Xv, (A-4) 

where u^ L is the four- velocity; x 11 is a um t spacelike vector orthogonal to it' 1 ; p is the energy density; p r is the 
pressure in the direction of x^j & n d is the pressure on the two-space orthogonal to x M . Since Tq\ = 0, it follows 
that the fluid is at rest in then system of reference of ([7])- Therefore, we can always choose u^ 1 = Sq /^/goo and 
X» = 5?/VFmj- Thus, 

T ° = p, Tl = - Pr , and T% = T$ = -p ± . (A-5) 

Equation of state: Now, from (|A-l[> - i|A-3|) . we find the equation of state of the effective matter. It is 

p = f(N,a)(p r + 2 P± ), (A-6) 

with 

m*) = A f- N \- 3N ^ N - 2 l ( A-7) 

4e(l-jV) + 3Na(N-2)' V 7 

where <r is given by ([5J) . Note that / = 1 for A" = and N — 2 the equation of state becomes 

p = Pr + 2p ± . (A-8) 

For JV = lwe get 

p = -(p r + 2 P± ). (A-9) 

In general, for any value of N the equation of state approximates (|A-8j) . for "small" values of a, i.e., for esil. 

Gravitational mass: In 4D, the gravitational mass inside a 3D volume V3 is given by the Tolman-Whittaker 
formula viz., 

M(r) = J (T ° - Tl - Tl - Ti)^g~ A dV 3 . (A-10) 
In isotropic coordinates, with goo = and gn — — e A(r \ this simplifies to 

M(r) = ir 2 e^+ A )/V, (A-ll) 

where the prime ' denotes derivative with respect to r. Using this equation, for the metric |(7J) we obtain 

2e I ar — 1 N 



M « = 7UI (A " 12) 



Here, we require e > in order to ensure the positivity of M . 
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